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We evaluate the ability of future cosmic microwave background (CMB) experiments to measure 
the power spectrum of large scale structure using quadratic estimators of the weak lensing deflection 
field. We calculate the sensitivity of upcoming CMB experiments such as BICEP, QUaD, BRAIN, 
ClOVER and Planck to the non-zero total neutrino mass indicated by current neutrino os- 
cillation data. We find that these experiments greatly benefit from lensing extraction techniques, 
improving their one-sigma sensitivity to by a factor of order four. The combination of data 
from Planck and the SAMPAN mini-satellite project would lead to a{M^) ~ 0.1 eV, while a value 
as small as (t(M^) ~ 0.035 eV is within the reach of a space mission based on bolometers with a 
passively cooled 3-4 m aperture telescope, representative of the most ambitious projects currently 
under investigation. We show that our results are robust not only considering possible difficulties 
in subtracting astrophysical foregrounds from the primary CMB signal but also when the minimal 
cosmological model (A Mixed Dark Matter) is generalized in order to include a possible scalar tilt 
running, a constant equation of state parameter for the dark energy and/or extra relativistic degrees 
of freedom. 

PACS numbers; 14.60.Pq, 95.35. +d, 98.80.Es 



o 



X 



I. INTRODUCTION 

Nowadays there exist compelling evidences for flavor 
neutrino oscillations from a variety of experimental data, 
that includes measurements of solar, atmospheric, reac- 
tor and accelerator neutrinos (for recent reviews, see e.g. 
P, 01) • The existence of flavor change implies that the 
three neutrinos mix and have non-zero masses, but os- 
cillation experiments only fix the differences of squared 
neutrino masses Am|]^ and Arn^i, which correspond to 
the values relevant for atmospheric (2.4 x 10"'^ eV^) and 
solar (8 x 10~^ eV^) neutrinos, respectively. 

Non-zero neutrino masses imply that the Cosmic Neu- 
trino Background (CNB), the sea of relic neutrinos that 
fill the Universe with a number density comparable to 
that of photons, influences the cosmological evolution in 
a more complicated way than that of a pure relativistic 
component. In particular, the contribution of the CNB 
to the present energy density of the Universe, measured 
in units of its critical value, is 



Pc 



93.14 /i2 eV 



(1) 



where h is the present value of the Hubble parameter in 
units of 100 km Mpc~^ and = mi + TO2 + is 
the total neutrino mass. From the experimental values of 



their mass differences, at least two neutrino mass states 
are non-relativistic today since both (Amii)i/2 _ 0.05 
eV and (Ato2i)^^^ 0.009 eV are larger than the present 
neutrino temperature T^, ~ 1.96 K ~ 1.7 x 10"'* eV. Since 
the current upper bound on Mj, from tritium decay ex- 
periments is of the order 6 eV (95% CL), we know that 
the neutrinos account for at least 0.5(1)% and at most 
50% of the total dark matter density, where the lower 
limit corresponds to the minimum of Mi, for masses or- 
dered according to a normal (inverted) hierarchy, charac- 
terized by the sign of Ato|]^. Thus, although in the first 
limit the cosmological effect of neutrino masses would be 
quite small, the minimal cosmological scenario is in fact 
a A Mixed Dark Matter (AMDM) model rather than a 
plain A Cold Dark Matter one. 



Considerable efforts are devoted to the determination 
of the absolute neutrino mass scale, which, combined 
with oscillation data, would fix the value of the light- 
est neutrino mass. The future tritium decay experiment 
KATRIN 4,j is expected to reach a discovery potential 
for 0.3 — 0.35 eV individual masses, while more stringent 
bounds exist from experiments searching for neutrinoless 
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double beta decay-'^. These will be improved in the near 
future iQ], but unfortunately they depend on the details 
of the neutrino mixing matrix. The quest for M^, will 
greatly benefit from cosmological observations, which of- 
fer the advantage of being independent of the neutrino 
mixing parameters since all flavors were equally popu- 
lated in the early Universe. 

Cosmology is sensitive to the neutrino masses through 
essentially two effects. First, the shape of the two-point 
correlation function -or power spectrum- of the Cosmic 
Microwave Background (CMB) temperature and polar- 
ization anisotropics on the one hand, and of the Large 
Scale Structure (LSS) mass density on the other, are 
both highly sensitive to the abundance of the various 
cosmological backgrounds: photons, baryons, cold dark 
matter, etc. The CNB is very specific in the sense that 
it behaves like a coUisionless relativistic medium at the 
time of acoustic oscillations before photon decoupling (at 
redshifts z > 1000), but like a non- relativistic fluid dur- 
ing most of structure formation (at redshifts z < 100, at 
least for one of the three neutrino mass states). There- 
fore, the CNB affects at least one of the three following 
quantities: the redshift of equality between matter and 
radiation; the redshift of equality between matter and 
dark energy; or the spatial curvature of the Universe. 
This effect can be observed in the CMB and LSS power 
spectra and its amplitude is at most of the order of (2/^) 
per cent (/^ is the current fraction of dark matter 
density in the form of neutrinos), which corresponds to 
only 1% in the limit ~ 0.05 eV. 

Fortunately, neutrino masses produce a second effect 
which is typically four times larger: on small scales neu- 
trinos do not cluster gravitationally because of their large 
velocities. Even today, the typical neutrino velocity of a 
non-relativistic eigenstate with mass is as large as 
V ~ 150 (1 eV/m^) km s~^. This simple kinematic effect, 
called neutrino free- streaming, is extremely important 
for the growth of non-relativistic matter perturbations 
(CDM and baryons) after photon decoupling. Indeed, 
the perturbation growth rate is controlled by the bal- 
ance between gravitational clustering and the Universe 
expansion. On small scales, free-streaming neutrinos con- 
tribute to the total background density p, but not to the 
total perturbation Sp, which shifts the balance in favor 
of the Universe expansion, leading to a smaller growth 
rate for CDM and baryon perturbations. This effect is 
of order (8/j/) per cent in the small-scale matter power 
spectrum 0, IE 111- 

There are various ways to measure the LSS power spec- 
trum. For instance, the galaxy-galaxy correlation func- 
tion can be obtained from galaxy redshift surveys, and 
the density perturbations in hot intergalactic gas clouds 



A claim of a positive signal exists |3 , which would correspond to 
an effective neutrino mass of order 0.1 — 0.9 eV. If confirmed, it 
would have a profound impact on cosmology. 



at redshift z 2 can be inferred from the Lyman-a forest 
region in the spectrum of distant quasars. At present, a 
total neutrino mass of 0.4 — 2 eV is disfavored at 95% CL 
[IS 111 111 111 III III El 111 111 E 113, 113, depending 
on the used CMB, LSS and other cosmological data. 

However, the most promising idea on the long term 
is to study the weak lensing effects induced by neigh- 
boring galaxy clusters. A lensing map can be recon- 
structed from a statistical analysis, based either on the 
ellipticity of remote galaxies or on the non-gaussianity of 
the CMB temperature and polarization anisotropy maps. 
Weak lensing offers several advantages. Unlike galaxy 
redshift surveys, it traces directly the total density per- 
turbation and does not involve any light-to-mass bias. 
Unlike Lyman-a forests data, it probes a large range of 
scales, which is particularly convenient for observing the 
step-like suppression of density perturbations induced by 
neutrino masses. In addition, weak lensing is sensitive 
to high redshifts, for which non-linear corrections appear 
only at very small scales. Finally, it enables tomographic 
reconstruction: by selecting the redshift of the sources, 
it is possible to obtain independent measurements of the 
power spectrum at various redshifts, in order to follow the 
non-trivial evolution of the spectrum amplitude caused 
by neutrino masses and/or by a possible evolution of the 
dark energy density. The best lever arm and the highest 
redshifts are encoded in the lensing of CMB maps, where 
the source is the photon last scattering surface, located at 
z ~ 1100, and the observed CMB patterns are sensitive 
to lenses as far as z 3 jH IH |2J| . In addition, CMB 
lensing observations do not require a devoted experiment: 
future CMB experiments designed for precision measure- 
ments of the primary CMB anisotropics offer for free an 
opportunity to extract lensing information. 

The first paper estimating the sensitivity of future 
cosmological experiments to small neutrino masses was 
based on the measurement of the galaxy-galaxy correla- 
tion function an analysis that was updated in Refs. 
25, 26, 27] and more recently in Ref. ,28]. The idea that 
weak lensing observations (from galaxy ellipticity) were 
particularly useful for measuring the neutrino mass was 
initially proposed in Ref. |2l|. Then, the first analysis 
based on CMB lensing extraction was performed in Ref. 
|30j , showing that an extremely small one-sigma error on 
the total neutrino mass -of the order of a{M,y) ~ 0.04 
eV- was conceivable for a full-sky experiment with a res- 
olution of 1 arc-minute and a sensitivity per pixel of 1 fiK 
for temperature, 1.4 /iK for polarization (these numbers 
were inspired from preliminary studies for the CMBpol 
satellite project). Soon after, Ref. jsl] studied the neu- 
trino mass sensitivity of future tomographic reconstruc- 
tions using, on the one hand, galaxy ellipticities in vari- 
ous redshift bins, and on the other CMB lensing, where 
CMB plays the role of the last redshift bin. The authors 
found that for sufficiently large cosmic shear surveys, it 
would not be impossible to reach ct(M^) ~ 0.02 eV. 

In this paper we want to come back to the prospects 
coming from CMB lensing alone, and try to improve the 
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pioneering analysis in 113,1^ in several directions. First, 
we analyze the potential of several CMB experiments ex- 
pected to produce results in the coming years, based on a 
realistic description of instrumental sensitivities. Second, 
we discuss the robustness of our results by analyzing (i) 
the consequences of simplifying assumptions in the con- 
struction of the Fisher matrix, (ii) the dependence of the 
final results on the accuracy of the foreground subtrac- 
tion process, and (iii) the impact of parameter degenera- 
cies which can appear when non-minimal cosmological 
scenarios are introduced. Finally, we study the sensitiv- 
ity of CMB experiments to the way in which the total 
neutrino mass is split among the three species. 



II. BASIC PRINCIPLES OF CMB LENSING 
EXTRACTION 

Weak lensing induces a deflection field d, i.e. a map- 
ping between the direction of a given point on the last 
scattering surface and the direction in which we observe 
it. At leading order ^| this deflection field can be writ- 
ten as the gradient of a lensing potential, d = V0. The 
(curl-free) deflection map and the lensing potential map 
can both be expanded in harmonic space 

0(n) = Y.<l>rYr{n) , (2) 
{dg±id^)(h) = ±z^d™±il^"(n) , (3) 

Im 

where n = (0, (f>) is a direction in the sky. There is a sim- 
ple relation between the deflection and lensing potential 
multipoles 

dr = -^^/l{^+T)cj,^ , (4) 

so that the power spectra Cf'^ = and C/"^ = 

are related through 

Cf'^ = /(/ + l)cf* . (5) 

In standard inflationary cosmology, the unlcnscd 
anisotropics obey Gaussian statistics in excellent approx- 
imation [3^ . and their two-dimensional Fourier modes 
are fully described by the power spectra C^^ where a and 
b belong to the {T, E,B} basis. Weak lensing correlates 
the lensed multipoles [22[ I3 according to 

{aTbf )CMB = {-irsl'5Z'cf + J2C{a,b)Z7'j!' cf,^' 

LM 

(6) 

where the average holds over different realizations (or 
different Hubble patches) of a given cosmological model 
with fixed primordial spectrum and background evolu- 
tion (i.e. fixed cosmological parameters). In this average, 
the lensing potential is also kept fixed by convention, 
which makes sense because the CMB anisotropics and 



LSS that we observe in our past light-cone are statisti- 
cally independent, at least as long as we neglect the inte- 
grated Sachs- Wolfe effect. In the above equation, Cf^ is 
the lensed power spectrum (which is nearly equal to the 
unlensed one, excepted for the B-mode power spectrum 
which is dominated, at least on small scales, by the con- 
version of E-patterns into B-patterns caused by lensing). 
The coefficients C(a, b)™iTj^' are complicated linear com- 
binations of the unlensed power spectra (7°^, C*;"" and 
C•f^ given in 

The qua dratic estimator method of Hu & Okamoto 
[35I l3a. [33 is a way to extract the deflection field map 
from the observed temperature and polarization maps. 
It amounts essentially in inverting Eq. ©. This is not 
the only way to proceed: Hirata & Seljak proposed an it- 
erative estimator method 3^ which was shown to be op- 
timal, but as long as CMB experiments will make noise- 
dominated measurements of the B-mode, i.e. at least for 
the next decade, the two methods are known to be equiv- 
alent in terms of precision. Even for the most precise ex- 
perimental project discussed in this work, the quadratic 
estimator method would remain nearly optimal (the last 
project listed in Table corresponds roughly to the hy- 
pothetical experiment called "C" in Ref. 

By definition, the quadratic estimator d{a, b) is built 
from a pair (a, b) of observed temperature or polarization 
modes, and its multipoles are given by the quadratic form 

d(a,6)f =AAf y^MTfi' aTbf , (7) 

where the normalization factor A/'£'' is defined in such way 
that d(a, b) is an unbiased estimator of the deflection field 

(d(a, 6)f )cMB - v/I(7TT)0f , (8) 

and the weighting coefficients W(a, b)Y^^iM minimize the 
variance of c?(o, b)'^ (which inevitably exceeds the power 
spectrum Cjf that we want to measure), i.e. minimize 
the coefficients a = a', 5 = 6' of the covariance matrix 

(d(a,6)f d(a',6')r)cMB = (-1)*'^5|;'< [Cf -|-iVf '''"'I . 

(9) 

Here the extra term Nf^"" * , which can be considered 
as noise, derives from the connected and non-connected 
pieces of the four-point correlation function {aba'b'). In 
Ref. js^, Okamoto & Hu derive a prescription for the 
weighting coefficients W(a, b)™J^iF such that the contri- 
bution of the connected piece is minimal, while that from 
the non-connected piece is negligible in first approxima- 
tion [33|. The weighting coefficients are rational func- 
tions of the observed power spectra Cf^ , Cf^ and Cf^, 
which include contributions from primary anisotropics, 
lensing and experimental noise. Therefore, if we assume 
a theoretical model and some instrumental characteris- 
tics, we can readily estimate the noise Nf""'^' expected 
for a future experiment. 

This method works for a given estimator d(a, b)^ un- 
der the condition that for at least one of the three power 
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spectra (C°^, C°", Cf*), the lensing contribution is much 
smaher than the primary anisotropy contribution. This 
is not the case for the pair ab = BB. Therefore, one 
can only build five estimators, for the remaining pairs 
ah e {TT, EE, TE, TB, EB}. The question of which one 
is the most precise heavily depends on the experimental 
characteristics. In addition, it is always possible to build 
a minimum variance estimator, i.e. an optimal combi- 
nation of the five estimators weighted according to the 



five noise terms Ni ° of the experiment under consid- 



jaha b 

eration. For the minimum variance estimator, the noise 
reads 



E 

.aba'b' 



aba b 
I 



-1 



(10) 



III. 



FORECASTING ERRORS WITH THE 
FISHER MATRIX 



For a future experiment with known specifications, it 
is possible to assume a cosmological fiducial model that 
will fit best the future data, and then to construct the 
probability L{x; 6) of the data x given the parameters 9 
of the theoretical model. The error associated with each 
parameter 9i can be derived from the Fisher matrix 



(11) 



computed in the vicinity of the best-fit model. Indeed, 
after marginalization over all other free parameters, the 
one-sigma error (68% confidence limit) on a parameter 
9i would be greater or equal to 



(12) 



In most cases, the forecasted errors depend only mildly on 
the exact values of fiducial model parameters; however, 
they can vary significantly with the number of free pa- 
rameters to be marginalized out, since complicated fidu- 
cial models with many physical ingredients are more af- 
fected by parameter degeneracies. 

It is usually assumed that for a CMB experiment cov- 
ering a fraction /sky of the full sky, the probability L of 
the data {a™} is gaussian, with variance C;. If the exper- 
iment observes only one mode, for instance temperature, 
then C; is just a number, equal to the sum of the fiducial 
model primordial spectrum and of the instrumental noise 
power spectrum. If instead several modes are observed, 
for instance temperature, E and B polarization, then C; 
is a matrix. Neglecting the lensing effect, we would get 



c, = 







CF 



NF 







NF 



. (13) 

where the C^^^s represent the power spectra of primary 
anisotropics (we recall that for parity reasons CF = 



CF 



0), and the N^^'s are the noise power spec- 



tra, which are diagonal because the noise contributing to 
one mode is statistically independent of that in another 
mode. It can be shown with some algebra that for any 
gaussian probability L, the Fisher matrix reads |39| 



1 f)C (9C 

= :jE(2^ + l)/^kyTrace[C'i — ] 



d9i 



(14) 



In fact, due to the lensing effect, the data is not ex- 
actly gaussian. However, the difference between the un- 
lensed and lensed power spectra for {TT, TE, EE) is 
so small that Eq. (|14|l remains approximately correct, 
at least when the B-mode is not included in the covari- 
ance matrix of Eq. H13|l . Beyond this issue, lensing of- 
fers the possibility to include an extra piece of informa- 
tion: namely, the map of the lensing potential -or equiv- 
alently, of the deflection vector- as obtained from e.g. 
the quadratic estimators method. Ideally, after lensing 
extraction, one would obtain four gaussian independent 
variables: the delensed temperature and anisotropy mul- 
tipoles (f™, i?™, B^), and the lensing multipoles dj". 
In this paper, we will take a fiducial model with no sig- 
nificant amplitude of primordial gravitational waves. In 
this case, the delensed B-mode is just noise and can be 
omitted from the Fisher matrix computation. Therefore 
the data covariance matrix reads 



Ci = 



CF 



CF 
CF 



NF 



CF 



CF 



nF 



cF 



(jdd ^ j^dd 



(15) 

where Cf'^ is the lensing power spectrum, Nf'^ the noise 
associated to the lensing extraction method (in our case, 
the minimum variance quadratic estimator) , and Cj'^ the 
cross-correlation between the unlensed temperature map 
and lensing map. This term does not vanish because of 
the late integrated Sachs- Wolfe effect: the temperature 
includes some information on the same neighboring clus- 
ter distribution as the lensing. Both Cf"* and Cj'^ can 
be computed numerically for a given theoretical model 
using a public Boltzmann code like CAMB 40], and then 
ISSf"^ can be calculated using the procedure of Ref. ^|. 
This computation can be performed in the full sky: in 
this work, we will never employ the flat-sky approxima- 
tion. Note that the B-mode does not appear explicitly 
in Eq. (|15|l . but actually information from the observed 
B-mode is employed in the two estimators (i(T, i?) and 
d{E,B). 

Using Eqs. H14|l and (|15|l . it is possible to compute a 
Fisher matrix and to forecast the error on each cosmo- 
logical parameter. Let us discuss the robustness of this 
method. There are obviously two caveats which could 
lead to underestimating the errors. 

First, we assumed in Eq. H15|) that the temperature and 
polarization maps could be delensed in a perfect way. In- 
stead, the delensing process would necessarily leave some 
residuals, in the form of extra power and correlations in 
the covariance matrix. However, this is not a relevant 
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issue, because we are using only the temperature and 
i?-polarization modes, for which the lensing corrections 
are very smalh therefore, considering a smaU residual or 
no residual at all makes no difference in practice. We 
checked this explicitly in a simple way. For a g iven the- 
oretical model, Boltzmann codes like CAMB [40| are able 
to compute both the lensed and unlensed power spectra. 
If the delensing process is totally inefhcient, we can say 
that unlensed temperature and polarization multipoles 
are recovered with an error of variance 

77TT \r<TT ATTi j^EE ^EEi 

El ^\Ci -Ci I, El =\Ci -Ci I, (16) 

that we can treat as additional noise and sum up to the 
Nf^ and Nj^^ terms in the matrix (|15|l . We checked 
numerically that even with such a pessimistic assump- 
tion, the final result does not change significantly, which 
is not a surprise since Ef"' <^ C^"'. We conclude that the 
assumption of perfect delensing performed in Eqs. H14|) 
and H15|l is not a problem in practice^. 

Second, we assumed a perfect cleaning of all the astro- 
physical foregrounds which contribute to the raw CMB 
observations. It is true that CMB experiments are op- 
erating in various frequency bands, precisely in order to 
subtract the foregrounds which frequency dependence is 
usually non-planckian. However, we still have a poor 
knowledge of many foregrounds, and some of them could 
reveal very difficult to remove, introducing extra non- 
gau ssianity and spoiling the lensing extraction process 
particular, the question of foreground sub- 
traction is related to the maximum / at which we should 
stop the sum in the Fisher matrix expression, i.e. to the 
smallest angular scale on which we expect to measure 
primary temperature and polarization anisotropics. If 
we assume a perfect cleaning, this value should be de- 
duced from instrumental noise. Beyond some multipoles 
(^max; ^max)j ^hc uoisc tcrms {Nf^, N^^) bccomc expo- 
nentially large. Thus, in practice, the sum in Eq. 1)14(1 
can be stopped at any / larger than both i^^^^ ^^'^ ^max- 
However, some foregrounds are expected to be impossi- 
ble to subtract on very small angular scales (e.g., radio 
sources, dusty galaxies, or polarized synchrotron radia- 
tion and dust emission), so for experiments with a very 
small instrumental noise, the covariance matrix could be 
dominated by foreground residuals at smaller I values 
than those where the instrumental noise explodes. 

Since we do not have precise enough data at high galac- 
tic latitude and on relevant frequencies, it is difficult at 



^ Note that replacing Cf^ by Cf^ in 1151 would actually be a mis- 
take. Indeed, in this case, the Fisher matrix would include the 
derivatives of the lensed power spectra with respect to the cos- 
mological parameter. So, the physical effect of each cosmological 
parameter on lensing distortions would be counted several times, 
not only in dCf'^/ddi but also in dCf/dOi, with a, 6 G {T, E}. 
This would introduce correlations which would not be taken into 
account self-consistently, and the forecasted errors would be ar- 
tificially small, as noticed in l30l . 



the moment to estimate how problematic foreground con- 
tamination will be, but it is clear that one should adopt 
a very careful attitude when quoting forecasted errors for 
future experiments with an excellent angular resolution. 
In the next sections, for each experiment and model, we 
will derive two results: one optimistic forecast, assum- 
ing perfect foreground cleaning up to the scale where the 
instrumental noise explodes (or in the case of the most 
precise experiments, up to the limit Z^axj ^max < 2750 
beyond which it is obvious that foreground contamina- 
tion will dominate); and one very conservative forecast, 
assuming no foreground cleaning at all. In that case, we 
take the foreground spectra Fj^'^ , Ff^ and Ff^ of the 
"mid-model" of Ref . j4j| , computed with the public code 
provided by the authors'^. This model is not completely 
up-to-date, since it is based on the best data available 
at the time of publication, and does not include impor- 
tant updates like the level of polarized galactic dust ob- 
served by Archeops on large angular scales J4]. Also, for 
simplicity, it assumes statistically isotropic and Gaussian 
foregrounds, with no TB or EB correlations. However 
this approach is expected to provide the correct orders of 
magnitude, which is sufficient for our purpose. We add 
these new terms to the covariance matrix of Eq. (|15|l , as 
if they were extra noise power spectra for the TT, EE 
and TE pairs. We consistently recompute Nf'^^ still us- 
ing the equations in Ref. |3^J but with these extra noise 
terms included, in order to model the worse possible loss 
of precision induce by foregrounds in the lensing extrac- 
tion process. We expect that the true error-bar for each 
cosmological parameter will be somewhere between our 
two optimistic and conservative forecasts. 



IV. EXPERIMENTAL SENSITIVITIES 

We consider seven CMB experiments which are rep- 
resentative of the experimental efforts scheduled for the 
next decade. The first two, based in the South Pole, 
are complementary: BICEP^ (Background Imaging of 
Cosmic Extragalactic Polarization) [i^ is designed for 
large angular scales, while QUaD^ (QUest at DASI, the 
Degree Angular Scale Interferometer) for small an- 
gular scales. The second experiment, which is already 
collecting data, is composed of the QUEST (Q and U 
Extragalactic Sub-mm Telescope) instrument mounted 
on the structure of the DASI experiment. A second set 
of experiments is scheduled in Antarctica at the French- 
italian Concordia station and in the Atacama plateau in 
Chile: the BRAIN^ (B-modes Radiation measurement 



For each experiment, we compute the foreground for each fre- 
quency channel, and then compute the minimum variance com- 
bination of all components. 

* http:/ /www. astro. caltech.edu/~lgg/bicep_front. html 
http: / / www. astro, cf. ac.uk/ groups / instrumentation /projects / quad / 

^ http:/ /apc-p7.org/APC_CS/Expcrienccs/Brain/index.phtml 
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from Antarctica with a bolometric INterferometer) p7l | 
instrument for measming large scales, and the ClOVER'' 
(CI ObserVER) instrument for intermediate scales. 
BRAIN and ClOVER are designed for unprecedented 
precision measurements of the S-mode for / < 1000. 
ClOVER was approved for funding by PPARC in late 
2004 and could be operational by 2008. At that time, 
the Planck^ satellite 49J should be collecting data: 
Planck has already been built and should be launched 
in 2007 by the European Space Agency (ESA). Beyond 
Planck, at least two space projects are under investiga- 
tion: the mini-satellite SAMPAN (SAtcUitc to Measure 
the Polarized ANisotropies) ^50,] for CNES (Centre Na- 
tional d'Etudes Spatiales) , and the more ambitious Infla- 
tion Probe project for NASA (National Aeronautics and 
Space Administration) , whose characteristics are not yet 
settled. The calculations of Ref. |3(| were based on num- 
bers inspired from preliminary studies for the CMBpol 
satellite project: a resolution of 1 arc-minute and a sen- 
sitivity per pixel of 1 fiK for temperature, 1.4 /iK for 
polarization. Here, the experiment that we will call In- 
flation Probe is based on one over many possibilities |5lj : 
a bolometer array with a passively cooled 3-4 m aperture 
telescope, with four years of multifrequency observations 
and a sensitivity of 2 /iK s^^/^ per channel. 

We list the expected instrumental performances of each 
experiment in Table Each instrument includes many 
detectors grouped in frequency bands or channels. In 
each channel, the detectors have a given spatial resolution 
described by the FWHM (Full- Width at Half-Maximum) 
Ob of the beam. For a given channel, one can estimate 
the temperature and polarization sensitivities per pixel of 
the combined detectors, and = A^- The channel 
noise power spectrum reads 



Nl"^: = {etAafexp + l)0bV81n2] 



(17) 



with a G {T,E,B}. The noise from individual channels 
can be combined into the global noise of the experiment 



Nr = 



aa\ — 1 



(18) 



Given this input, the computation of the lensing noise 
N^"^ can be performed numerically following Ref. [35l| . 
In Fig. n we show our results for the lensing noise of 
each experiment, based on each quadratic estimator and 
on the combined minimum variance estimator. In Fig. [21 
we gather information on the noise for the TT, EE and 
dd power spectra for each experiment. The error-bars 
AC;"" displayed in Fig. |21 include both cosmic variance 
and instrumental noise, and assume a multipole binning 



http:/ /www-astro. physics. ox. ac.uk/~act/clover. html 

* http:/ /sci . esa. int / science-e / www / area/index. cfm? fareaid=17] 
aiid ,http:/ /www.planck.fr/. 
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0.29 


0.41 
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0.42 


0.59 


150 
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220 
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TABLE I: Sensitivity parameters of the CMB projects 


con- 


sidered in this work: /sky is the observed fraction of the sky, 


V the center frequency of the channels 


m GHz, 


the FWHM 


(Full- Width at Half-Maximum) in arc-minutes 


, At the tem- 


perature sensitivity per pixel in /iK and A_b 


= Ab the po- 


larization sensitivity. For all experiments, we 


assumed one 


year of observations, except for the Inflation Probe sensitiv- 


ity based on four years. 








of width A^ = 7 until ^ 70, and then A/ ^ 






^^'" = \/(2/ + 1)'a^/.,.(^'" + ^''')- 


(19) 



The top graphs in Figs. & |21 correspond to the 
BICEP-fQUaD and BRAIN-hClOVER combinations. 
Computing the Fisher matrix for each pair of experiments 
is not a trivial task, due to the different sky coverages. 
We follow a method which is certainly not optimal, but 
has the merit of simplicity. Since in each case, one ex- 
periment is optimized for large scales and the other for 
smaller scales, we assume that below a given value Ic all 
multipoles are evaluated from BICEP or BRAIN only, 
while for I > Ic they are taken from QUaD or ClOVER. 
In Eqs. I|14l I15|l . this amounts in considering /sky as a 
function of I, and in replacing /sky(Oi Nf'^ and Nf'^ 
by their BICEP/BRAIN value for I < l^, or by their 
QUaD/ClOVER value for / > 1^. The lensing noise Nf"^ 
is then computed for the combined experiment, following 
the same prescriptions. For each pair of experiments, we 
optimized the value of Ic numerically by minimizing the 



7 




10"- 



10""^ 



10"' 



10"= 



10"" 









TT 


Plancl^^ 


EE 




TE 




TB 




_EB 




mv 







FIG. 1: For six CMB experiments or combinations of experiments, we show the expected noise power spectrum Ni for the 
quadratic estimators d{a, b) built out of pairs ab £ {TT, EE, TE, TB, EB}, and for the combined minimum variance estimator 
(mv). The thick line shows for comparison the signal power spectrum Cf'' — The sum of the two curves Nf'^ + Cf'* 

represents the expected variance of a single multipole ci(a,6)[". 



forecasted error on the total neutrino mass Mi,. In both 
cases, we found that I ^ 300 was optimal. This method 
might be less favorable for BRAIN+CIOVER than for 
BICEP+QUaD, because the first pair of experiments has 
a large overlap in /-space, for vifhich multipoles could be 
derived from the t-wo combined datasets. 



impressively small noise power spectrum Nf'^. QUaD has 
both an excellent resolution and a very good sensitivity, 
and should provide an extremely precise measurement of 
T and E modes on small angular scales. Therefore, the 
three quadratic estimators d{T,T), d{E,E) and d{T,E) 
are particularly efficient. 



We find that BICEP+QUaD is able to reconstruct the 
lensing multipoles in the range 2 < I < 200 with an 



The main goal of the BRAIN+CIOVER combined ex- 
periment is to improve the determination of the i?-mode 
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FIG. 2: For the same six CMB experiments or combinations of experiments as in figure|21 we show the expected binned error 
on the reconstructed power spectra: from top to bottom, Cf'^ (using the minimum variance quadratic estimator), Cf^ and 
Cf^. The curves represent the power spectra of the fiducial model described in section IVI 



performed by BICEP+QUaD, especially on large and in- 
termediate scales {I < 1000), which are particularly im- 
portant for detecting gravity waves from inflation. This 
should be achieved with a sensitivity which is even bet- 
ter than that of BICEP and QUaD, but at the expense 
of a poorer resolution in the case of ClOVER, leading 
to large errors for small-scale polarization. In total, this 
design is roughly equivalent to that of BICEP+QUaD 
in terms of lensing extraction: BRAIN+CIOVER is also 



able to reconstruct the lensing multipoles in the range 
2 < Z < 200. The best estimator is now d{E,B), known 
to be particularly useful, since E and B are correlated 
only due to lensing. In this sense, future lensing deter- 
minations by BRAIN-f ClOVER and by BICEP+QUaD 
can be seen as complementary, and therefore both par- 
ticularly interesting. 

The Planck satellite has a resolution comparable to 
QUaD, but a poorer sensitivity than the last four exper- 
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iments. This explains why the lensing noise shown in 
Fig. n looks a bit disappointing: the signal marginally 
exceeds the noise only around ^ ~ 40. However, we 
should keep in mind that Planck will observe the full 
sky (which leads to /sky — 0.65, once the galactic cut 
has been taken into account), while BICEP+QUaD or 
BRAIN+CIOVER explore only smaU regions. There- 
fore, for a given I, Planck makes many more indepen- 
dent measurements of multipoles (T™, EJ^), and conse- 
quently, also of d™. In Fig. [21 one can check that Planck 
still makes a more precise determination of the lensing 
power spectrum than BICEP-|-QUaD: both experiments 
are able to constrain Cf"^ up to Z ^ 1100, but the satellite 
provides smaller errors. 

Since Planck is not very sensitive to i?-modes, and 
BRAIN is limited by its small sky coverage, there will 
be room after these two projects for improving i3-mode 
observations on large angular scales, in view of observing 
inflationary gravitational waves. This would be the tar- 
get of the SAMPAN mini-satellite project, which would 
be a full-sky experiment with excellent sensitivity but 
poor resolution. We find that for the minimum vari- 
ance estimator, the noise Nf'^ would be at the same level 
for Planck and SAMPAN. However, it is interesting to 
note that Sampan has a good d{E, B) estimator, while 
Planck is better with d{T, T). Therefore, it sounds par- 
ticularly appealing to combine the two full-sky experi- 
ments, that is technically equivalent to assuming a super- 
experiment with twelve channels (nine from Planck and 
three from SAMPAN). The resuhs (in the fifth graphs 
of Figs, n] & 0) show that with such a combination one 
could lower the noise Nf'^ by a factor two for the mini- 
mum variance estimator, in order to constrain Cf"^ up to 
I - 1300. 

Finally, the (hypothetical) version of the Infiation 
Probe satellite that we consider here has an extremely 
ambitious resolution and sensitivity, such that the instru- 
mental error would be better than cosmic variance for the 
_B-mode until I ~ 1500. For such a precise experiment, 
assumptions concerning foreground subtraction play a 
crucial role, since it is very likely that foreground resid- 
uals will start dominating the observed power spectrum 
before instrumental noise. The last graphs in Figs.n&El 
which assume perfect foreground cleaning up to Z 2500, 
show that lensing multipoles could be recovered up 
to to / 800, while Cf'^ could be constrained up to at 
least I - 2500. 



V. FUTURE SENSITIVITIES TO NEUTRINO 
MASSES 

For each experiment, we compute the Fisher matrix 
following Eqs. (^El, for a AMDM fiducial model with 
the parameter values as given below, and considering 
two possibilities for the number of free parameters that 
should be marginalized out. 

The first possibility is the minimal alternative on 




-1 - 



10 100 1000 

I 

FIG. 3: Logarithmic derivatives of the lensing power spec- 
trum Cf with respect to each cosmological parameter. The 
derivatives with respect to ojb and tJm have been rescaled in 
order to fit inside the figure. 



the basis of current observations: we marginalize 
over eight free parameters, standing for the current 
baryon density c^b = f^b^^j the current total mat- 
ter density Wm — f^m , the current dark energy 
density SI a, the total neutrino mass M^, in eV, the 
primordial curvature power spectrum amplitude 
and tilt ti^, the optical depth to reionization t and 
the primordial helium fraction yuc, to which we as- 
sign the values (wb, Wm, f^A, M^, ln[10^°A5], Us, r, ?/hc) = 
(0.023,0.143,0.70,0.1,3.2,0.96,0.11,0.24). We assume 
no spatial curvature and tensor contribution. Note 
that the reduced Hubble parameter derives from h = 

The second possibility, describing non-minimal phys- 
ical assumptions, is to marginalize over three extra pa- 
rameters: the scalar tilt running a, which can be non- 
negligible in some inflationary models with extreme as- 
sumptions; the dark energy equation-of-state parame- 
ter w] and finally, extra relativistic degrees of free- 
dom which would enhance the total radiation density, 
parametrized by the effective number of neutrino species 
Ncff (for instance, A'eff = 4 means that the Universe con- 
tains a background of extra relativistic particles with the 
same density as one extra massless neutrino species). 
In the fiducial model, these parameters take the val- 
ues {a,w, Ncs) = (0,-1,3). Our purpose is to find out 
whether such extra free parameters open up degeneracy 
directions in parameter space, that would worsen the sen- 
sitivity to neutrino masses. It has been shown in recent 
analyses that these parameter degeneracies indeed ap- 
pear with current CMB and LSS data (see 14, 15] for 
iVoff and mm for w). 

The derivative of the lensing power spectrum with 
respect to each of these eleven parameters are shown in 
Fig. 121 with the exception of the derivatives with respect 
to r and ?/ho which are null, and with respect to Ag which 
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Free parameters: 




8 parameters 


of minimal AMDM 




same + {a, 






Lensing extraction: 


no 


no 


J CO 


J CO 


no 


no 


yes 


yes 


Foreground cleaning: 


perfect 


none 


p6rf6ct 


I10I16 


perfect 


none 


perfect 


none 


QUaD+BICEP 


1.3 


1.6 


0.31 


0.36 


1.5 


1.9 


0.36 


0.40 


BRAIN+CIOVER 


1.5 


1.8 


0.34 


0.43 


1.7 


2.0 


0.42 


0.51 


Planck 


0.45 


0.49 


0.13 


0.14 


0.51 


0.56 


0.15 


0.15 


SAMPAN 


0.34 


0.40 


0.10 


0.17 


0.37 


0.44 


0.12 


0.18 


Planck+SAMPAN 


0.32 


0.36 


0.08 


0.10 


0.34 


0.40 


0.10 


0.12 


Inflation Probe 


0.14 


0.16 


0.032 


0.036 


0.25 


0.26 


0.035 


0.039 



TABLE II: Expected 1-a error on the total neutrino mass Af^ in eV for various CMB experiments or combinations of them. 
The first (last) four columns correspond to a AMDM model with eight (eleven) free parameters. For each of the two models, the 
four columns show the cases with or without lensing extraction, and with two extreme assumptions concerning the foreground 
treatment: perfect subtraction or no subtraction at all. 



is trivial. All derivatives were computed using the public 
Boltzmann code CAMB enabling the highest accu- 
racy options and increasing the accuracy_boost parame- 
ter to five. Whenever possible, we evaluated double-sided 
derivatives, and searched for optimal step sizes such that 
the results were not affected by numerical errors (from 
the limited precision of the code) nor by contributions 
from higher-order derivatives. 

We quote the results for the total neutrino mass in 
Table im assuming either eight or eleven free parameters. 
For each of the two cases, we compare the forecasted er- 
rors with and without lensing extraction, i.e. using either 
a2x2ora3x3 data covariance matrix, in order to eval- 
uate the impact of the extraction technique. Finally, in 
each of the four sub-cases, we quote the results obtained 
assuming perfect foreground cleaning or no cleaning at 
all, in order to be sure to bracket the true error. Should 
we trust more the results based on the eight or eleven 
parameter model? This depends on future results from 
cosmological observations: in absence of strong observa- 
tional motivation for extra parameters, one will probably 
prefer to stick to the simplest paradigm; however, the 
next years might bring some surprises, like for instance 
the detection of a variation in the dark energy density. 

Let us comment the results for each experiment. The 
combination QUaD-|-BICEP benefits a lot from lensing 
extraction, since the error decreases from approximately 
1.5 eV to at least 0.4 eV. These results are found to be 
robust against foreground residuals and extra parameter 
degeneracies. It is interesting that with QUaD+BICEP 
it should soon be possible to reach in a near future - 
using CMB only- the same precision that we have today 
combining many observations of different types (galaxy- 
galaxy correlation function, Lyman-a forests) which are 
affected by various systematics. The situation is al- 
most the same for BRAIN-I-CIOVER, which should also 
achieve a{M,y) ~ 0.4 eV using lensing extraction. 

Planck should make a decisive improvement, lower- 
ing the error to cr(M^) ~ 0.15 eV, in excellent agree- 
ment with the results of Ref. [s^. Note that without 
lensing extraction the error would be multiplied by three 
(by four in the case with extra free parameters). We do 



not find a significant difference between the forecasted 
errors in the eight and eleven parameter models. SAM- 
PAN alone is slightly more efficient than Planck, and 
the combination Planck-|-SAMPAN is the first one to 
reach a{M,j) ~ 0.1 — 0.12 eV, even in the pessimistic case 
of large foreground residuals and extra free parameters. 
Thus these future CMB lensing data could help in break- 
ing the parameter degeneracy between Afj, and w j54| . 
that would still be problematic at the level of precision 
of Planck (without lensing extraction) combined with 
the galaxy-galaxy correlation function extracted from the 
Sloan Digital Sky Survey. 

Finally, the version of the Inflation Probe satellite that 
we consider here is able to reach <j{M^) — 0.035 eV both 
in the eight and eleven parameter cases. Note that when 
we take instead the CMBpol specifications of Ref. [30l| . 
we exactly reproduce their forecast a{M^) = 0.04 eV (de- 
rived for an intermediate case with ten parameters). It 
is interesting to see that even with such a precise experi- 
ment, the results are robust against foreground contami- 
nation, since in absence of any cleaning the forecast error 
increases only by 15%. 

We show in Fig.^jthe correlation between and each 
free parameter of the eleven-dimensional model, in the 
cases of Planck and Inflation Probe, with and without 
lensing extraction. In the parameter basis used in this 
work, the neutrino mass appears to be mainly degenerate 
with ujmi and the lensing extraction process removes most 
of this degeneracy. 



VI. FUTURE SENSITIVITIES TO THE 
NEUTRINO MASS SPLITTING 

In principle, the LSS power spectrum is not sensitive 
only to the total mass M^, but also to the way in which 
the mass is distributed among the three neutrino states. 
The reason is twofold: the amount by which the gravi- 
tational collapse of matter perturbation is slowed down 
by neutrinos on small scales depends on the time of the 
non-relativistic transition for each eigenstate, i.e. on the 
individual masses; and the characteristic scales at which 
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FIG. 4: 1-a confidence Hmits on the pairs (M^, 6i), for each 
parameter 9i in our eleven-dimensional model. The red solid 
(green dashed) contours are those expected for Planck (Infla- 
tion Probe). For each case, the smaller (larger) ellipse corre- 
sponds to the forecasts with (without) lensing extraction. 



significant mass, or two degenerate states in the Inverted 
Hierarchy case (IH). In general, for a given one can 
calculate the difference between the matter power spec- 
trum in the two cases, as has been computed numerically 
in Ref. ^J. 

We would like to study whether the lensing power spec- 
trum derived from a very precise CMB experiment like 
Inflation Probe would be able to discriminate between 
the two models. For this purpose, we take the eight pa- 
rameter model of section fVl and complete it with a ninth 
parameter: the number of massive neutrinos jv™'^'''"™, 
which could be equal to 1, 2 or 3 (the remaining species 
are assumed to be exactly massless). In a NH scenario 
with Ml, > 0.1 eV, the mass of the third neutrino is 
not completely negligible: so, we expect the difference 
between our simplified scenario with ^^^^^^^'^ — i and 
that with TV"^^*'^'™ = 2 to be more pronounced than the 
difference between realistic NH and IH scenarios (assum- 
ing the same total mass A/^ in all models). This state- 
ment is confirmed by the numerical results of Ref. [28| . 
So, if we could show that an experiment like Inflation 
Probe will be unable to differentiate between the sketchy 

^massive ^ ^ ^^^j ^massive ^ 3 modcls, WC WOuld COn- 

clude that a fortiori it will not discriminate between the 
NH and IH scenarios. 

We repeated the computations of section with 
a ninth free parameter jv™'*'''"™ with fiducial value 
^massive ^ jsj^^g ^Yia.t the parameter TV^nassivo g^ould 

not be confused with the total effective neutrino num- 
ber iVoff , which was a free parameter in the last section, 
and remains fixed to iVcff = 3 in the present one. We 
found for Inflation Probe - including lensing extraction 
and assuming perfect foreground cleaning- a one-sigma 
error a (tv™''""^™) = 2.8. We conclude that the experi- 
ments and techniques discussed in the present paper are 
far from sufficient for discriminating between the NH and 
IH scenarios. In any case, as shown in Ref. future re- 
sults on the total neutrino mass from very precise cosmo- 
logical data should be interpreted in a slightly different 
way for the NH and IH cases. 



VII. CONCLUSIONS 



the free-streaming effect of each neutrino family is im- 
printed in the power spectrum depends on the value of 
the wavelengths crossing the Hubble radius at the time 
of each non-relativistic transition, i.e. again on the indi- 
vidual masses. 

The neutrino masses are differently distributed among 
the three states in the two possible mass schemes, or hi- 
erarchies, as shown e.g. in Fig. 1 of "2^. For a total mass 
Ml, larger than 0.2 — 0.3 eV all neutrino states approxi- 
mately share the same mass mo = M^/3, in the so-called 
degenerate region. Instead, for smaller the splitting 
between the individual masses is more visible, and for 
the minimum values of one finds that in the Normal 
Hierarchy case (NH) there is only one neutrino state with 



We have studied the ability of future CMB experiments 
to measure the power spectrum of large scale structure, 
using some quadratic estimators of the weak lensing de- 
flection field. We inferred the sensitivity of these ex- 
periments to the non-zero neutrino masses indicated by 
neutrino oscillation data. Our aim was to extend the 
pioneering paper by Kaplinghat, Knox & Song [s^l by 
further investigating several directions. 

First, we based our analysis on the following list of 
forthcoming CMB experiments (either operational, ap- 
proved or still in project): BICEP, QUaD, BRAIN, 
ClOVER and Planck, SAMPAN and Inflation Probe, 
taking into account their detailed characteristics. We 
found that even before Planck, ground-based experi- 
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ments should succeed in extracting the lensing map with 
good precision, and could then significantly improve the 
bounds on neutrino masses. We also found that the 
SAMPAN mini-satellite project would be able to reduce 
the Planck error a{M^) from approximately 0.15 eV to 
0.10 eV. Finally, the hypothetical version of Inflation 
Probe that we considered would reach a spectacular sen- 
sitivity of a{M^) ~ 0.035 eV. 

We also tried to discuss two questions raised by the 
analysis of Ref. [s^]: first, is it really accurate to base the 
Fisher matrix computation on perfectly delensed maps on 
the one hand, and on the reconstructed lensing map on 
the other? Second, is it realistic to estimate the noise 
variance of the lensing quadratic estimators without tak- 
ing into account any residual foreground contamination? 
Our answer to these two questions is positive: we did not 
provide an exact treatment of these very technical issues, 
but we tried to systematically bracket the results be- 
tween two over-optimistic and over-pessimistic assump- 
tions, and concluded that the error forecast method of 
Ref. is robust. 

Finally, we investigated the issue of parameter degen- 
eracies involving the neutrino mass, by comparing the 
results in a simpler model than that of Ref. [s^l with 
those in a more complicated one. Our extended cos- 
mological model allows for a scalar tilt running, a dark 
energy equation of state parameter w ^ —1, and extra 
degrees of freedom parametrized by the effective number 
of massless neutrinos iVeff- These extra parameters were 
not chosen randomly. The tilt running was shown to be 
slightly degenerate with the neutrino mass in an analy- 
sis involving current CMB and LSS data The same 
holds for the equation of state of dark ener gv [541 and for 
the effective number of massless neutrinos |l4l |l5| . How- 
ever, our results indicate that future CMB experiments 
will be able to resolve these degeneracies, since we do not 
find significant discrepancies between the neutrino mass 
errors obtained for our two cosmological models. 

Fortunately, CMB lensing extraction should be re- 
garded as only one of the most promising tools for mea- 



suring the absolute neutrino mass with cosmology. It 
could be combined with future data from tomographic 
galaxy cosmic shear surveys, which will be very sensitive 
to neutrino masses 'si^ . The cross-correlation of LSS in- 
formation with CMB temperature anisotropics could also 
reveal very useful for the purpose of measuring M^, |55j . 
In the method employed in the present paper, the cor- 
relation between temperature and lensing (the Td term) 
is already taken into account, but it affects the final re- 
sults only marginally. More interesting should be the 
cross-correlation of future data from large cosmic shear 
surveys with that from CMB anisotropies. 

In conclusion, our results show that there are good 
perspectives to detect non-zero neutrino masses using 
future CMB lensing data, since even in the less favor- 
able case of the smallest ~ 0.05 eV in the NH mass 
scheme the Inflation Probe experiment alone could make 
a marginal detection (between the one and two sigma 
levels). Obviously the sensitivity is enhanced for larger 
values of il/^, in particular for the mass degenerate and 
quasi-degenerate regions but also for the minimum of 

~ 0.1 eV in the IH case. The information on 
from analyses of cosmological data will be complemen- 
tary (and vice versa) to the efforts in terrestrial projects 
such as tritium beta decay and neutrinoless double beta 
decay experiments. Of course any positive result on the 
absolute neutrino mass scale will be a very important in- 
put for theoretical models of particle physics beyond the 
Standard Model. 



Acknowledgments 

We would like to thank Karim Benabed, Simon Prunet 
and Jonathan Rocher for extremely useful discussions. 
This work was supported by a MEC-IN2P3 agreement. 
SP was supported by the Spanish grants BFM2002-00345 
and GV/05/017 of Generalitat Valenciana, as well as by 
a Ramon y Cajal contract of MEC. 



[1] M. Maltoni, T. Schwetz, M.A. Tortola and J.W.F. Valle, 

New J. Phys. 6, 122 (2004) hep-ph/0405172 . 
[2] G.L. Fogli, E. Lisi, A. Marrone and A. Palazzo, Prog. 

Part. Nucl. Phys., in press hep-ph/0506083 . 
[3] K. Eitel, Nucl. Phys. Proc. Suppl. 143, 197 (2005). 
[4] A. Osipowicz et al. [KATRIN Coll.], hep-ex/0109033 
[5] H.V. Klapdor-Kleingrothaus et al, Phys. Lett. B 586, 

198 (2004) hep-ph/0404088 . 
[6] S.R. Elliott and P. Vogel, Ann. Rev. Nucl. Part. Sci. 52, 

115 (2002) hep-ph/0202264 . 
[7] J. Lesgourgues and S. Pastor, in preparation. 
[8] J.R. Bond, G. Efstathiou and J. Silk, Phys. Rev. Lett. 

45, 1980 (1980). 
[9] W. Hu, D.J. Eisenstein and M. Tegmark, Phys. Rev. 

Lett. 80, 5255 (1998) astro-ph/9712057 . 
[10] D.N. Spergel et al [WMAP Coll.], Astrophys. J. Suppl. 



148, 175 (2003) astro-ph/0302209 . 
[11] S. Hannestad, JCAP 0305, 004 (2003) 

astro-ph/0303076 . 
[12] M. Tegmark et al, Phys. Rev. D 69, 103501 (2004) 

astro-ph/0310723 . 
[13] V. Barger, D. Marfatia and A. Tregre, Phys. Lett. B 595, 

55 (2004) hep-ph/0312065 . 
[14] S. Hannestad and G. Raffelt, JCAP 0404, 008 (2004) 

hep-ph/0312154 . 
[15] P. Grotty, J. Lesgourgues and S. Pastor, Phys. Rev. D 

69, 123007 (2004) hep-ph/0402049 . 
[16] R. Rebolo et al, Mon. Not. Roy. Astron. Soc. 353, 747 

(2004) astro-ph/0402466 . 
[17] U. Seljak et al, Phys. Rev. D 71, 103515 (2005) 

astro- ph/0407372 . 
[18] G.L. Fogli et al, Phys. Rev. D 70, 113003 (2004) 



13 



hep-ph/0408045 . 
[19] K. Ichikawa, M. Fukugita and M. Kawasaki, Phys. Rev. [38 
D 71, 043001 (2005) astro-ph/0409768 . 

C. J. MacTavish et al, astro-ph/0507503 [39 
A.G. Sanchez et al, Mon. Not. Roy. Astron. Soc. 366, 
189 (2006) astro-pli/0507583 [40 
F. Bernardeau, Astron. Astrophys. 324, 15 (1997) 
astro-ph/9611012 . [41 
M. Zaldarriaga and U. Seljak, Phys. Rev. D 58, 023003 

(1998) astro-ph/9803150 . [42 
U. Seljak and M. Zaldarriaga, Phys. Rev. Lett. 82, 2636 

(1999) astro-ph/9810092 . [43 

D. J. Eisenstein, W. Hu and M. Tegmark, Astrophys. J. 
518, 2 (1998) astro-ph/9807130 . [44 
J. Lesgourgues, S. Pastor and S. Prunet, Phys. Rev. D 
62, 023001 (2000) hep-ph/9912363 . [45 
S. Hannestad, Phys. Rev. D 67, 085017 (2003) 
astro-ph/02 11106. [46 
J. Lesgourgues, S. Pastor and L. Perotto, Phys. Rev. D 
70, 045016 (2004) hep-ph/0403296 . [47 
A.R. Cooray, Astron. Astrophys. 348, 31 (1999) 
astro-ph/9904246 . [48 
M. Kaplinghat, L. Knox and Y.S. Song, Phys. Rev. Lett. 
91, 241301 (2003) astro-ph/0303344 . [49 
Y.S. Song and L. Knox, Phys. Rev. D 70, 063510 (2004) [50^ 
astro-ph/0312175 . 

CM. Hirata and U. Seljak, Phys. Rev. D 68, 083002 [51 
(2003) astro-ph/0306354^. [52 
N. Bartolo, E. Komatsu, S . Matarrese and A . Riotto, 
Phys. Rep. 402, 103 (2004) 'astro-ph/0406398 . [53 
U. Seljak, Astrophys. J. 463, 1 (1996) 
astro-ph/9505109 . 

T. Okamoto and W. Hu, Phys. Rev. D 67, 083002 (2003) [54 
astro-ph/0301031 . 

W. Hu, Astrophys. J. 557, L79 (2001) [55 
astro-ph/0 105424 . 

W. Hu and T. Okamoto, Astrophys. J. 574, 566 (2002) 



astro-ph/0111606 . 

A. Cooray and M. Kesden, New Astron. 8, 231 (2003) 
astro-ph/0204068 . 

M. Tegmark, A. Taylor and A. Heavens, Astrophys. J. 
480, 22 (1997) astro-ph/9603021 . 

A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. 
538, 473 (2000) astro-ph/9911177 . 

A. Amblard, C. Vale and M.J. White, New Astron. 9, 
687 (2004) astro- ph/0403075 . 

A. Cooray, M. Kamionkowski and R.R. Cald well, Phys. 
Rev. D 71, 123527 (2005) astro-ph/0503002' . 

M. Tegmark, D.J. Eisenstein, W. Hu and A. de Oliveira- 
Costa, Astrophys. J. 530, 133 (2000) astro-ph/9905257 . 
N. Ponthieu et al, Astron. Astrophys. 444, 327 (2005) 
astro-ph/0501427 . 

B. C. Keating et al, in Polarimetry in Astronomy, ed. by 
S. Fineschi, Proceedings of the SPIE 4843, 284 (2003). 
S. Church et al [QUaD CoU.], New Astron. Rev. 47, 1083 

(2003) . 

M. Piat et al [BRAIN-CLOVER Coll.], SF2A 2004 pro- 
ceedings (2004) astro-ph/0412590 . 

B. Maffei et al. [BRAIN-CIOVER CoU.], EAS Publica- 
tions Series 14, 251 (2005). 
J.A. Tauber, Adv. Space Res. 34, 491 (2004). 
F.R. Bouchet et al, SF2A 2005 proceedings (2005) 
astro-ph/0510423 . 
J. Bock, private communication. 

M. Bowden et al, Mon. Not. Roy. Astron. Soc. 349, 321 

(2004) astro-ph/0309610 . 

The Planck Consortia, "The Scientific Programme of 
Planck 2004 (Blue Book)", ESA Pubhcation Division 

(2005) . 

S. Hannestad, Phys. Rev. Lett. 95, 221301 (2005) 
astro-ph/0505551 . 

K. Ichikawa and T. Takahashi, |astro-ph/0510849| 



